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Abstract — We study the characteristics and radiation mecha- 
nism of antenna superstrates based on closely located periodical 
grids of loaded wires. An explicit analytical method based on 
the local field approach is used to study the reflection and 
transmission properties of such superstrates. It is shown that 
as a result of proper impedance loading there exists a rather 
wide frequency band over which currents induced to the grids 
cancel each other, leading to a wide transmission maximum. 
In this regime radiation is produced by the magnetic dipole 
moments created by circulating out-of-phase currents flowing 
in the grids. An impedance matrix representation is derived for 
the superstrates, and the analytical results are validated using 
full-wave simulations. As a practical application example we 
study numerically the radiation characteristics of dipole antennas 
illuminating finite-size superstrates. 

Index Terms — Superstrate, wire grid, capacitively loaded wire, 
local field, gain enhancement 

I. Introduction 

Artificial dielectrics, implemented e.g. as arrays of wires, 
have been used for decades as light-weight beam shaping 
elements, see e.g. [l]-[6] for some early contributions and 
[7], [8] for some more recent results. It is nowadays well 
understood that some electromagnetic (photonic) crystals sim- 
ulate the behavior of homogeneous materials with ultra-low 
refractive index near stop-band edges [9], [10], [11]. This 
observation has lead to several studies concerning directive 
antenna design using such periodical structures to shape the 
beam of a low-gain primary radiator, e.g. [12]— [14]. It is also 
known that a slight local change in the period of electromag- 
netic crystals leads to localized resonant modes which can 
be used for the realization of devices radiating energy in a 
very narrow angular range [15]. Typically, radiating sources 
are placed inside Fabry-Perot resonant cavities formed by 
removal of rows of wires or dielectric rods in a periodical 
lattice [16]— [21]. The long -known ideas of creating resonant 
cavities using partially reflecting sheets [22] have also been 
widely employed in recent gain enhancement designs [23]- 
[27]. Similarly, frequency-selective surfaces on top of ground 
planes have been used to implement Fabry-Perot resonant 
cavities, e.g. [28]. Antenna superstrates are another well- 
known gain enhancement tool [29], [30]. As in the case of 
defected electromagnetic crystal structures [16]— [21] one of 
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the key ideas of different superstates is to allow radiation from 
a primary source to spread over a larger radiating aperture. 
Recently, very exotic structures such as material covers or 
superstrates aimed to possess double-negative behavior [31] 
have been proposed for this purpose [32], [33], [34]. 

The goal of this work is to present an analytical methodol- 
ogy for studying the characteristics of a class of composite 
structures proposed as superstrates for low-gain antennas. 
More precisely, we consider two closely located parallel 
grids of capacitively loaded wires having possibly a grid of 
continuous (purely inductive) wires in between them. The 
methodology presented here can, however, be extended to 
study the characteristics of other kinds of superstrates based 
on multiple grids of wires loaded with arbitrary distributed 
impedance. The physical mechanism responsible for the oper- 
ation of the proposed superstrates is clarified. We also briefly 
discuss the extension of the methodology to handle strongly 
coupled dipole grids. The transmission properties of such 
grids have in the past been analyzed with an application 
to frequency-selective surfaces or electromagnetic band-gap 
structures, e.g., in [36], [37], [38]. In these works the response 
has been calculated using an iterative method based on solving 
a set of coupled integral equations for the unknown element 
currents and aperture fields. The method presented here is 
explicit, and despite taking into account all the near field 
interactions between the grids (higher-order Floquet modes) it 
remains computationally efficient and physically illustrative. 

The novelty of the proposed superstrate designs is to excite 
resonant oscillations in capacitively loaded wires (this way 
creating a wider radiating aperture) while maintaining partial 
transparency of the structure. In the designs considered here 
the source field induces out-of-phase currents to the loaded 
wire pairs, and the pairing currents give rise to resonant 
magnetic dipole moments that produce aperture radiation. 
After heuristically explaining the operational principle of the 
proposed superstrates we move on calculating the reflection 
and transmission properties of the superstrates. The analytical 
results are validated using full-wave simulations, and for a 
physical illustration we calculate the electric and magnetic 
dipole moments induced to one of the structures. It is shown 
that strong resonant enhancement of the superstrate aperture 
field can be achieved with very thin structures. To further 
clarify the physics behind the superstrate operation we derive 
an impedance matrix representation directly from the electro- 
dynamic model. Finally, a high-gain antenna structure utilizing 
the proposed superstrates is implemented in a commercial 
method-of-moments solver. 
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II. Heuristic operational principle 

The superstate studied in this work is schematically de- 
picted in Fig. QJ. The first and third grids (see the numbering 
of grids in Fig. QJ) consist of parallel capacitively loadeqj 
wires, while the second grid consists of continuous (inductive) 
thin wires. Please note that in Fig. Q] round wires have 
been replaced by strips to reflect one feasible implementation 
technique. However, in the following analytical derivation we 
assume round wires since electrically thin wires and strips 
behave in a similar manner. To account for the difference in 
the geometry an effective radius of wires is introduced in the 
analytical analysis. 

The properly designed structure is expected to show an 
analogy with double-negative behavior: In a certain frequency 
range the metal strips constituting a capacitively loaded wire 
can be interpreted as short dipoles, and every pair of such 
dipoles (sitting in grids 1 and 3) shows a magnetic resonant 
mode which can be utilized to create a response similar to 
that of split rings [39], [40]. On the other hand, the long 
conductors provide plasma-like electric response as in wire 
media [41], [42]. Recently, it has been reported that only short 
wire pairs are enough to obtain double-negative behavior [43], 
[44]. For describing the radiation properties of the superstate 
we concentrate, however, on the distributed dipole moments 
induced in the structure. 

Consider first only a pair of two capacitively loaded wires 
(assuming that grid 2 is absent). At a certain frequency the 
currents induced to individual pairing metal strips cancel each 
other allowing the incident wave to propagate through the 
grids. To widen the pass-band one can bring a grid of solid 
wires in between the loaded wire grids. In this case two 
modes appear due to the capacitances and mutual coupling 
between the loaded wires and the continuous wires. In other 

'in this work capacitive loading is created by cutting the wires. 
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Fig. 1 . a) A schematic illustration of the proposed superstrate (not necessarily 
in scale), b) Unit cell geometry. 



words, there are now two possibilities for the currents in 
the grids to cancel each other, and by properly choosing the 
unit cell dimensions both transmission bands can be merged 
together. Since the superstrate is volumetric (though the total 
thickness 2h is in practise extremely small compared with 
the wavelength) the cancelation of total electric current does 
not mean the cancelation of radiation: The amplitude of 
the reflection coefficient is minimized at this frequency (and 
consequently the amplitude of the transmission coefficient is 
maximized), however, the magnetic dipole moments created 
by the pairing out-of-phase currents produce radiation. 

Another important practical issue, discussed below in more 
detail, is the interaction between a closely located primary 
radiator and the loaded wire grids. It will be shown that in case 
of capacitively loaded wires the source can be brought very 
close to the grids without deteriorating the radiation resistance. 
This allows us to design low-profile and effectively radiating 
structures. 

III. Local field method to study the superstrate 

PROPERTIES 

A. Plane-wave reflection and transmission coefficients and 
induced electric and magnetic dipole moments 

At this point we assume that the grids in the superstrate are 
infinite in the z and x directions. Moreover, the transversal 
thickness t of the inclusions (equivalently the radius tq of 
round wires) is small compared to the wavelength. We assume 
the normal plane-wave incidence (though the analysis can be 
carried out for oblique incidence as well) and write the external 
electric field as 



E c 



u z E cxt e~ jky . 



(1) 



We write down the equation for the local electric field acting 
on the surface of a reference element in each of the grids. 
Consider first the unloaded reference wire, and assume that 
the wires are thin enough so that we can represent them as 
infinitely long current lines located at the wire's axes. The 
relation between the amplitude of the local electric field E l ° c 
acting on the wire surface and current / flowing on the wire 
surface reads 

(2) 



where ceo is the susceptibility of the wire which can be 
determined from the boundary condition on the wire surface 
(see e.g. [45]): 

-l 



a 



(3) 



where rj is the wave impedance in the matrix material (host 

(2) 

material for the wires), and H Q ' represents the Hankel func- 
tion of the second kind and zero order. If the wires are loaded 
with a certain uniformly distributed impedance per unit length 
Z, the inverse of the wire susceptibility becomes [35], [46]: 



a 



Z. 



(4) 



The amplitude of the local electric field acting on the reference 
wires can also be expressed in the following manner: 

aT l h = E c 



/3(0)/i + P(h)I 2 + 0(2h)h, 



(5) 
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E ext e- jkh + (3(h)h + (3(0)I 2 + (3(h)h 



(6) 



E c 



-> 2kh + (3(2h)h + (3(h)I 2 + /3(0)/ 3 , (7) 



where Jj (i = 1, 2, 3) is the amplitude of current induced to 
the grid i, and /3(0) is so called self interaction coefficient 
(eq. (11) in [41], see also [47]): 
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and (3(h), f3(2h) are mutual interaction coefficients (eqs. (12) 
and (13) in [41]): 
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Above 7 w 0.5772 is the Euler constant, and is the x- 
component of the n-th Floquet mode wave vector [41]: 

*4 n) - -jV(Wd) 2 ~ fc2 > MVOI > 0. (10) 

Physically (3(0) tells how strongly the wires located in the 
same grid as the reference wire influence the local field, 
whereas (3(h), (3(2h) measure the influence of the other two 
grids. Now we can cast eqs. (|5])-(|7| in the following form: 
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From the above equation currents in different grids can be 
solved by simple inversion. The averaged current density 
induced to the grids is defined as 



7 _ h 

J% — — 7 5 * — 1,2, 3. 
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(12) 



In the far zone, the reflected field of one grid is a plane- 
wave field, and in the grid plane this plane-wave field can be 
expressed as [46] 

/; - 

(13) 



Ef c = —lj t . 
i 2 

After this we can calculate the total reflection and transmis- 
sion coefficients. The reflection and transmission coefficients 
referred to the plane of the first grid read: 

^ 3 1 £ ,sc e _: '' c ' l ( i_1 ) 
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As an example, we calculate the dipole moments induced 
to the structure consisting of two grids of capacitively loaded 
wires (in this case grid 2 is absent in Fig.QJ). The electric and 



magnetic dipole moments confined in the unit volume (Fig.QJ)) 
are calculated as [48] 



m 



up dV, 



fi / r x J dV, 



(16) 



(17) 



where p is the charge density in the unit volume, and J 
is the corresponding volumetric current density. When the 
integrations are conducted in our geometry the dipole moments 
take the following form: 



P : 



Id 



(h + h), 



hU , T , 



(18) 



(19) 



1 ) Discussion on the modeling approach: It is known that 
at low frequencies (frequencies at which the structural dimen- 
sions are small compared to the wavelength) the dispersion 
properties of lattices of capacitively loaded wires and dipole 
scatterers rather closely resemble each other [35]. Noticeable 
deviation is seen only close to the frequencies where the 
strip length in the loaded wires becomes resonant. At these 
frequencies (and at higher frequencies), for accurate analysis 
one should model the first and third grids as grids of dipole 
scatterers. The principal ideas of the above derivation would 
still be the same, however, different models for inclusion 
polarizabilities and interaction constants should be used [49]- 
[51]. Thus, a similar local field approach as introduced in the 
present paper can also be used to treat grids in which the 
unit inclusions can be modeled as small dipole scatterers. 
It will be demonstrated below that in the frequency regime 
where the superstrates proposed in this paper are feasible for 
antenna beam shapers, the capacitively loaded wire model 
leads to accurate results. It is also important to note that the 
method described above can be used to study the properties 
of superstrates formed by an arbitrary number of grids, and 
the distributed impedance loading the wires can be arbitrary. 

2) Numerical validation: For the validation of the ana- 
lytical method we choose for the superstate the structural 
parameters considered in [34] (the notation is clear from 
Fig. 0): I = 12.11 mm, l d = 10.11 mm, t = 1 mm, d = 4.8 
mm, h = 1 mm. We introduce equivalent round-cross-section 
wires, and use text book formulas to estimate the distributed 
capacitance due to the slits in the wires in first and third grids. 

For calculating the reflection and transmission coefficients 
we consider separately two structures. In the first one only 
the two grids of capacitively loaded wires are present. In the 
second case we bring the grid of continuous wires in between 
the loaded wire grids as shown in Fig. QJ. Reflection and 
transmission coefficients are calculated using eqs. (Tufl i and 
(TTSb . We also simulate the coefficients using Ansoft HFSS 
[52] (the unit cell for HFSS simulations is shown in Fig. QJ)). 
The results of calculations and simulations are depicted in 
Figs. |2] and [3] We observe that the agreement between the 
analytical calculations and simulations is good up to 13 GHz. 
At this frequency the length of the metal strips in the loaded 
wires is already roughly A/2, so it is clear that at frequencies 
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higher than 13 GHz (in this particular case) the local field 
model based on loaded wires fails. Nevertheless, for antenna 
applications we are interested in the frequency regime where 
the superstrates are partially transparent (frequency regime in 
which the reflection minima and transmission maxima lie, see 
Figs. |2] and [3]), and at these frequencies the analytical model 
works well. 

The amplitudes of electric and magnetic dipole moments 
induced in the unit volume of the superstrate consisting of 
two grids of loaded wires are shown in Fig. |4] As speculated 
earlier, the magnetic dipole moment shows a strong resonance 
in the vicinity of the frequency where the plane-wave reflection 
coefficient is minimized. It is worth to note that at the fre- 
quency of magnetic resonance the thickness of the superstrate 
is only A/12. 

B. Impedance matrix representation 

Based on the analysis of induced currents we derive an 
impedance matrix representation [53] connecting the total 
(averaged) electric field in the grid planes to the averaged 
current densities flowing in the grids: 





jgtot ^ 










V 


4 ot ) 





Z u 


Z\2 




fj 2 ) 




Z21 


Z22 


£) 




■ (20) 


Z31 


Z32 




v i ) 





Such representation is useful, e.g., when representing the grids 
as equivalent discrete loads in a waveguide modeling free 
space plane-wave propagation [3]. The expressions deduced 
below allow also to represent the mutual interaction between 
the grids at the impedance description level. In the quasi- 
static limit the deduced impedance expressions can be used 
to construct an equivalent circuit for the superstrates. For the 
sake of clarity we conduct below the derivation for the first 
grid in details. 

Using the impedance matrix notation the total field in the 
plane of the first grid reads 
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Fig. 3. Calculated and simulated reflection and transmission coefficients for 
two grids of capacitively loaded wires. 



Following the definition of the impedance matrix [53], Z\\ 
is the total electric field in the first grid plane divided by the 
averaged current flowing in that grid, assuming that currents in 
the other grids are zero. The component of the total (averaged) 
electric field parallel to the grids, defined in the plane of grid 
i is by definition [46] 



(22) 



When we solve E ext from (0 (now I2 = I3 = 0) and substi- 
tute the result in (1221 . the following expression is obtained: 



= -d 0(0) 
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Next, consider the definition of Z\2- The plane-wave term 
disappears in d22"b since now the current flowing in the first 
grid is by definition zero [53]. Thus, we get (after we solve 




Fig. 2. Calculated and simulated reflection and transmission coefficients for 
two grids of capacitively loaded wires with one grid of continuous wires in 
between them. 



Fig. 4. Dashed line: Electric dipole moment induced in the unit volume (in 
[Asm]), solid line: Magnetic dipole moment induced in the unit volume (in 
[Vsm]). The superstrate consists of two grids of capacitively loaded wires. 
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E ext from Q assuming that 1% = I3 = 0): 



Similarly for Z 13 : 



-d(3{h)J 2 = Z 12 J 2 . 



h=i 3 =o 



(24) 
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= -d(3(2h)J 3 = Z 13 J 3 . 

I 1= I 2 =0 

After repeating the above procedure for the second and third 
grid we get the following expressions for the impedances: 
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Z23 = -P(h)d, 
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IV. Finite size superstrates illuminated by wire 

ANTENNAS 

In this section we study and discuss the performance of 
finite-size superstrates when illuminated by practical antennas. 
Throughout the rest of the paper our primary source is a 
resonant A/2-dipole antenna. Naturally, the superstrates could 
also be illuminated with other types of antennas, such as with 
microstrip antennas, for example. 

A. Large superstrate with infinitely long wires 

Before designing an actual prototype antenna we study how 
well the key ideas for the superstrate performance, outlined 
with the plane-wave excitation, apply when the superstrates 
are illuminated with a resonant dipole, and the number of 
wires in the superstrates is finite. Namely, we study possible 
cancelation of electric currents in such structures. To avoid 
the standing wave interference caused by a finite wire height, 
and more clearly see the possible cancelation effect, we 
assume in this subsection that the wires are infinitely long. 
To ease the physical interpretation we consider here only 
the structure containing two grids of loaded wires, since in 
this case there is only one possibility (one frequency) for the 
cancelation of the total electric current, please see Fig. [3] Later, 
when designing the prototype antenna we consider also the 
superstrate containing a grid of unloaded wires. 

The number of wires in the grids is fixed to be 20, and 
we also fix the source position, see Fig. [5] The structural 
parameters are the same as introduced in Section III A. To 
handle the infinite wire length we use the method described 
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Fig. 5. Geometry of the structure used to study the distribution of current 
induced to pairing wires (not necessary in scale). 
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Fig. 6. Amplitudes of currents induced to the wires at different frequencies. 
Black dots denote the amplitude of current induced to wires located in the 
grid closest to the dipole (Fig. [5), white dots denote the amplitude of current 
induced to the wires in the other grid, a) / = 11.6 GHz, b) / = 10.4 GHz, 
c) / = 13.4 GHz. 



in [54] to solve the currents induced to the grids by taking 
into account the dipole excitation and the mutual interaction 
between all the wires. The distribution of induced current 
(amplitude and phase) is studied at three different frequencies, 
namely at 11.9, 10.4, and 13.4 GHz. Based on the result in 
Fig. [3] we predict that 11.9 GHz corresponds to the frequency 
at which the amplitude of current induced to the pairing wires 
is more or less the same, but the phase difference between the 
pairing currents is roughly ir. 

The amplitudes of the currents induced to the pairing wires 
are shown in Fig. [6] (the feed current amplitude equals unity). 
The phase differences of currents induced to the pairing wires 
are depicted in Fig. [7] We observe that at 11.9 GHz the 
amplitudes of currents induced to the pairing wires are close 
to the same value in most of the pairs. Slight difference is 
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Fig. 7. Phase difference (normalized to tv) between currents induced to the 
pairing wires at different frequencies, a) / = 11.6 GHz, b) / = 10.4 GHz, 
c) / = 13.4 GHz. 
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seen only in the wires closest to the dipole, but this difference 
is compensated by the fact that also the phase difference 
between the currents induced to the corresponding wire pairs 
slightly differs from tt. In those wire pairs where the current 
amplitudes are practically the same, also the phase difference 
is practically tt. The aforementioned implies that, indeed, the 
total electric current induced to the grids cancels as was 
speculated already in Section II. At 10.4 GHz the amplitudes 
of induced currents are much lower than at 1 1.9 GHz. Also the 
amplitudes of several pairing currents considerably differ from 
each other, and in many pairs the phase difference between the 
currents considerably differs from tt. Even though the phase 
differences of pairing currents is close to tt in all the pairing 
wires at 13.4 GHz (Fig. [7};), the total electric current is not 
canceled at this frequency since the amplitudes of currents 
constituting the pairs differ from each other in many pairs, 
Fig.B. 

B. High-gain antenna structure utilizing the proposed super- 
states 

In this subsection we design and implement a high-gain 
prototype antenna in a commercial method-of-moments solver 
FEKO [55]. The schematic illustration of the studied antenna 
structure is depicted in Fig. [5] We consider separately two 
structures: In the first case the grid of unloaded wires is absent 
(when analyzing the results we call this case the "2 grids" 
case), and in the second case the structure is as depicted in 
Fig-H("3 grids"). The radiation characteristics of both antenna 
structures are studied at the three frequencies introduced in the 
previous subsection. 

The dipole length is tuned to be A/2 at every considered 
frequency. The wire length is fixed to be l w = 3 A at 11.9 
GHz. Also the width (in ^-direction) of the superstrate is 
approximately 3A at 11.9 GHz. To enhance the directivity 
of the antenna structure, and make it more attractive as a 
possible high-gain solution, we create a resonant cavity using 
a metal reflector behind the grids, Fig. [8] Some of the available 
earlier contributions concerning the implementation of Fabry- 
Perot resonant cavities using partially reflecting sheets or 




other structures are listed in the Introduction. The distance 
between the metal reflector and the first wire grid equals 2d 
which roughly corresponds to A/2 at 11.9 GHz. The other 
relevant structural parameters are the same as introduced in 
connection with the plane-wave analysis (see Section III A)@ 
In the following the horizontal plane (H-plane) corresponds 
to the plane perpendicular to the wires, and vertical (E-plane) 
corresponds to the plane parallel to the wires. 

1) Example results: For this study the interesting antenna 
characteristics are the antenna radiation pattern, the input 
impedance Z- m , and the maximum boresight directivity £) max . 
We have not performed any systematic optimization of one or 
more of these characteristics, thus, from many points of view 
the presented antenna performance might not be optimal. The 
reader is also reminded that compared to the situation analyzed 
in Section III where the number of wires is infinite in the 
x-direction and the wire length is infinite, the finite number 
of wires and the finite wire height lead to the formation of 
standing waves along x and z-directions. This standing wave 
formation depends also on the location of the source, and it 
rather strongly influences the radiation characteristics. 

The simulated radiation patterns are depicted in Figs. |9TfTT1 
and the main simulated radiation characteristics are gathered in 
Table U Even without systematically optimizing the structural 
parameters we observe a good boresight directivity at 11.9 
GHz, which is consistent with the analysis presented so far 
in this paper. For comparison, without the superstrate the 
boresight directivity of a metal-plate backed dipole is 8.5 dB. 
In addition to this, the input resistance (in case of lossless 
structures, the radiation resistance) is high and makes the 
antenna easy to match to a standard 50 Ohm cable. It was ana- 
lytically shown in [56] that when a wire antenna is positioned 
in close vicinity of capacitively loaded wires, the currents 
induced to the closest wires add up constructively with the 

2 In FEKO we model round wires, and the capacitance caused by slits in 
the strips is modeled as an equivalent bulk capacitor. 
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Fig. 8. Schematic illustration of the prototype antenna geometry (not 
necessary in scale). The non-colored rods denote capacitively loaded wires, 
the grey rods denote unloaded wires. The black rod is the resonant A/2-dipole. 
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Fig. 9. Radiation patterns at 11.9 GHz. 
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TABLE I 

Simulated radiation characteristics. 
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Fig. 10. Radiation patterns at 10.4 GHz. 



source current, and therefore the radiation resistance remains 
high. This, indeed, appears to be the case. For comparison, in 
case when all the wires depicted in Fig. [8] are unloaded, the 
simulated input resistance at 11.9 GHz is only 0.7 Ohm. 

At 10.4 and 13.4 GHz the radiation characteristics are not 
so favorable for a feasible high-gain solution. At 10.4 GHz 
the peak directivity is weaker with both superstrates partly 
because the resonant cavity is not tuned for this particular 
frequency, and partly because the effective superstrate aperture 
is smaller than at 11.9 GHz. We see from Fig. |2]that the plane- 
wave reflection coefficient of the 3-grid superstrate is still 
very small at 10.4 GHz due to the dual resonance behavior. 
The corresponding effect is possibly seen as a much more 
narrow beamwidth of the main lobe (in boresight direction) 
as compared with the 2-grid structure. However, the finite 
structural dimensions clearly cause diffraction in the backward 
and sidelobe directions, and therefore the obtained maximum 
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directivity is only moderate. We believe, however, that after 
optimizing the structure at 10.4 GHz the directivity could be 
enhanced (if desired antenna performance would be required 
over a wider frequency band). At 10.4 GHz and also at 
13.4 GHz the input resistance is visibly smaller than at 
11.9 GHz, and this might lead to difficulties in achieving 
high radiation efficiencies when realistic losses are taken into 
account (especially with the 2-grid structure at 10.4 GHz). 
Additional simulations (not shown) on the current distribution 
induced to the wires in the 2-grid superstrate indicate that at 
11.9 GHz there is only one current maxima along the wires 
(in the z-direction, Fig. [8}, and the current amplitude in all the 
wires (excluding the two pairs of wires closest to the dipole) is 
rather close to the same value. This means that the superstrate 
aperture is rather uniformly excited, and a good boresight 
directivity is expected. At 10.4 GHz only every second pair is 
strongly excited. This is also the case at 13.4 GHz, in addition 
to this, at 13.4 GHz standing waves with three current maxima 
in z-direction are induced to the wires, and this is seen as the 
beam splitting in the E-plane. 

V. Conclusion 

We have presented a method to analyze the characteristics 
of a class of antenna superstrates based on periodical grids 
of loaded and unloaded wires. In addition to this, the physics 
behind the radiation mechanism of such superstrates has been 
clarified. The local field method has been used to solve the 
currents induced to the superstrates by the incident plane- 
wave field, and the reflection and transmission properties of 
the superstrates have been studied. The analytical results have 
been verified using numerical simulations. As a practical ap- 
plication example we have considered finite-size superstrates 
illuminated by dipole antennas, and discussed the effects of 
finite structural dimensions to the superstrate performance. A 
high-gain antenna structure utilizing the proposed superstrates 
has been implemented in a commercial full-wave simulator, 
and the results show promising antenna performance. 
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